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ABSTRACT 

We have imaged the central ~ 1 kpc of the circumnuclear starburst disk in the galaxy NGC 253 in 
the HCN(l->0), HCO+(1^0), and CO(1->0) transitions at 60 pc resolution using the Owens Valley 
Radio Observatory Millimeter- Wavelength Array (OVRO). We have also obtained Atacama Pathfinder 
Experiment (APEX) observations of the HCN(4^3) and the HCO+(4^3) lines of the starburst 
disk. We find that the emission from the HCN(1— >0) and HCO + (l— >0) transitions, both indicators 
of dense molecular gas, trace regions which are non-distinguishable within the uncertainties of our 
observations. Even though the continuum flux varies by more than a factor 10 across the starburst 
disk, the HCN/HCO + ratio is constant throughout the disk, and we derive an average ratio of 1.1±0.2. 
From an excitation analysis we find that all lines from both molecules are subthermally excited and 
that they are optically thick. This subthermal excitation implies that the observed HCN/HCO" 1 " line 
ratio is sensitive to the underlying chemistry. The constant line ratio thus implies that there are no 
strong abundance gradients across the starburst disk of NGC253. This finding may also explain the 
variations in £hcn/^hco+ between different star forming galaxies both nearby and at high redshifts. 
Subject headings: ISM: molecules - galaxies: individual (NGC253) - galaxies: starburst - galaxies: 
ISM 



1. INTRODUCTION 

Constraining the physical properties of the dense 
molecular gas phase is important for understanding the 
relation between the gas properties and the star for- 
mation process in nearby and high-z galaxies. HCO + 
has successfully been used as a tracer of dense gas 
(riH 2 > 10 4 cm~ 3 ) bot h in star forming regio ns of 
the "Milky Way (e.g., IChristopher et all 120051) and 
in nearby galaxies (e.g.. iNguven-Q-Rieu et alT 119921 : 
iBrouillet et al.1 [2005f) . Recently, it has also been de- 
tected at high redsh ifts towards the Cloverleaf quasar 
([Riechers et all 120061) and towards APM 08279+5255 
( Garcia-Burillo et al.ll2006l ). HCN is another good tracer 
of de nse gas and has also been observed both locally 
(e.g., iGao fc Solomon! I2004D and towards high redshift 
sourc e s~(e.g.. iSolqmon et al.l 120031 iVanden Bout et al.l 
2004) iCarilli et alj [20051 iWagg et alj [20051 IGao et al.l 



200l . IGao fc Solomon! ([20041 ) found that the HCN lu- 



minosity is correlated with the far-infrared l uminosity 
for sta rforming galaxies over several decades. IWu et all 
(2005) found that this correlation holds all the way to 
the smaller scales of the local star forming regions, i.e. 
over 7-8 orders of magnitude. 

The number of nearby galaxies in which molecular 
gas observations allow resolving the emission across the 
dense gas starburst regions is limit ed to only a few. A t 
a distance of only about 3.5 Mpc (jRekola et al.l l2005h . 
NGC 253, a member of the Sculptor group of galax- 
ies, is one of the closest starbursting galaxies. This 
makes it one of the best-suited laboratories for stud- 
ies of a circumnuclear starburst and its effect on the 
ambient interstellar medium (ISM). Some early detec- 

1 Max-Planck-Institut fur Astronomie, Konigstuhl 17, D-69117 
Heidelberg, Germany 

2 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, 
D-53121 Bonn, Germany 

3 Department of Astronomy and Radio Astronomy Laboratory, 
University of California, 601 Campbell Hall, Berkeley 



tions of H CN and HCO+ t oward s the center of NGC253 
came from iRickard et al.l ([1977[ ) and IRickard fc Palmer] 
( 19811) . both using single-dish telescopes. iPaglione et al.1 
(1995) presented first interferometry maps of HCN to- 
wards NGC253 using the Nobeyama Millimeter Array. 

We here present Owens Valley Radio Observatory 
(OVRO) interferometry maps and Atacama Pathfinder 
Experiment (APEX) 4 single-dish mapping of both HCN 
and HCO + , which we use for studying the spatial distri- 
bution and excitation of the two molecules in the circum- 
nuclear starburst of NGC253. The paper is structured 
as follows: section 2 presents the observations and data 
reduction, the results are given in §3, and in §4 we dis- 
cuss the implications both in the context of other nearby 
galaxies and for high redshift sources. At a distance of 
3.5 Mpc, 1" on the sky corresponds to 17.0 pc. 

2. OBSERVATIONS AND REDUCTION 

2.1. OVRO observations 

We observed NGC 253 in the HCN(l->0) (u lest = 
88. 63185GHz), HCO+(1^0) [v ieBt = 89.18852GHz) and 
CO(1->0) (f rost = 115.2712GHz) transitions with the 
OVRO Millimeter- Wavelength Array in its C and E con- 
figurations. The HCN and HCO + observations were ob- 
tained simultaneously during four tracks taken between 
1999 November and 2000 February (60 channels, with a 
channel width of 4MHz, corresponding to ~13.5kms _1 ); 
the CO observations were taken during six tracks from 
2001 October to 2002 March (120 channels with a width 
of 2MHz ~5.2kms _1 ); the duration of one track was 
typically 5 hr. The observations for each array were 
edited and calibrated separately with the MMA and the 
miriad packages. The uu-data were inspected and bad 
data points due to either poor atmospheric conditions 

4 This publication is based on data acquired with the Ata- 
cama Pathfinder Experiment (APEX). APEX is a collaboration 
between the Max-Planck-Institut fur Radioastronomie, the Euro- 
pean Southern Observatory, and the Onsala Space Observatory. 
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Fig. 1. — Channel maps of the HCN (thick contours) and HCO+ (thin contours) emission in NGC 253. The beamsize for both maps is 
identical (3.9" X 2.7") and two contours are plotted for each cubes at 50 and 200 mjybeam - 1 (la rms = 22 mjybeam -1 ). The central 
velocities for each channel map are given in the top left corner (in units of kms -1 ). 



or shadowing were removed, after which the data were 
calibrated. Flux calibration was determined by observ- 
ing Neptune during most observing runs. For observ- 
ing runs without planet observations the flux calibration 
has been bootstrapped using the gain calibrator 3C454.3 
(~9 Jy). These calibrators and an additional noise source 
were used to derive the complex bandpass corrections. 
The nearby calibrators B0135-247 (-0.37 Jy), B0048- 
097 (-0.72 Jy), and B0003-066 (-1.25 Jy) were used as 
secondary amplitude and phase calibrators. We estimate 
our flux calibration to be accurate to within 15%. 

For each transition, we produced a data cube using 
"robust" weighting which was cleaned down to a level of 
twice the rms noise. This resulted in beam sizes (and 
position angles) of 3.79"x2.56" (13?7) and 3.84"x2.64" 
(8?9) for the HCN and the HCO+ observations, re- 



spectively. As one of the main aims of this study is 
to compare the properties of the HCN and the HCO + 
emission, we convolved both cubes to a common res- 
olution of 3.9" x 2.7" (10?0), which resulted in an rms 
noise of — 22mJybeam~ 1 in one 4 MHz channel. The 
resulting beam size for the CO(l— ►()) observations is 
4.9"x2.7" (— 1?9) with an rms of — 65mJybeam _1 per 
2 MHz channel. The channels not containing HCN and 
HCO + line emission were combined to create a contin- 
uum map at 3 mm; the resulting rms in this image is 
5mJybeam . We note that the resolution of the CO 
and HCN maps i s comparable to that of th e interfero- 
metric maps from iPaglione et al.l (|1995l . l2004h , while the 
HCO + ha s a higher resolut ion than the Hat Creek Array 
map from lCarlstroml (|1990l ). Furthermore, the resolution 
of our 3 mm continuum map is a bit better than those of 
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iPeng et all {T996) and lGarci'a-Burillo et all (|2000l ). 

As the channel width of the CO(l— >0) observations is 
different from the HCN (l->0) and HCO+(1^0) obser- 
vations, we regridded the CO(l— >0) data cube to match 
the channel width of the other transitions. Subsequently, 
we used the high signal-to-noise ratio CO(l— >0) data 
cube to define, in each channel map, the regions that 
contain emission (this is called the "master" cube). We 
used this master cube to blank both the HCN and HCO+ 
cubes. All integrated maps shown here are based on these 
blanked cubes and contain the same areas of emission for 
all three transitions. 

2.2. APEX observations 

We also obtained observations of the HCN(4— >3) 
(i/ rest = 354.50547 GHz) and HCO+(4->3) (^ rcst = 
356.73413 GHz) transition using the APEX 12m tele- 
scope on Chajnantor. For these frequencies the half- 
power beamwidth of the antenna is 18". The observa- 
tions were carried out during two runs in 2006 September 
in good weather conditions (typical water columns: 0.7- 
lmm) using the dual-side-band receiver APEX-2A. Cal- 
ibration was achieved using the standard hot-cold-sky 
observing scheme roughly every 10 minutes. For both 
runs typical system temperatures were ~ 200 K (T*). 
Spectra were recorded using the Fast Fourie r Trans- 
form Spectrometer on APEX (|Klein et al.ll2006| ) provid- 
ing a bandwidth of 1 GHz with an effective resolution 
of 120kHz (0.12kms _1 ). The telescope pointing was 
checked on the nearby carbon star RScl and found to 
be accurate to ~ 3". The observations were performed 
in the standard raster-mapping mode covering the full 
region where emission has been detected in the OVRO 
maps in both transitions. Data processing was done us- 
ing the CLASS software package. In the data process- 
ing, we dropped all scans with distorted baselines, sub- 
tracted linear baselines from the remaining spectra, and 
then rebinned to a velocity resolution of 20kms _1 . Fi- 
nally, temperatures were converted to T m b using a for- 
ward efficiency an d a main beam effi ciency of 0.97 and 
0.73 respectively ([Glisten et al .1120061 ). The conversion 
from T m b to flux is 32 JyK -1 at 350 GHz. We estimate 
our flux calibration to be accurate to within 15%. 

3. RESULTS 

Figure Q] shows the OVRO channel maps for the 
HCN(l^O) and HCO+(1^0) data. It is striking that 
their spatial distributions are essentially identical. Both 
molecules are tracing dense molecular gas in the circum- 
nuclear starburst region of NGC 253. 

Figure [2] shows the integrated maps of HCN(l^O) 
and HCO + (l— >0) line emission, respectively. The 
slightly different morphology towards the southwest in 
the HCO+(1^0) map compared to HCN(l^O) is likely 
due to the slightly lower signal-to-noise ratio in the map. 
The HCN(4^3) and HCO+(4->3) spectra from APEX 
are shown in Figure [3] The line profiles are very similar 
at all observed positions, and the rotation of the circum- 
nuclear gas is nicely visible in both transitions. The in- 
tegrated fluxes for all transitions are given in Table Q] 
along with the line l uminosities, which we derive using 
the equations given in lSolomon et al.l (|1997| ). The values 
for HCN and HCO + grou nd transitions are consistent 
with previous results from iNguven-Q-Rieu et all (|1989h 
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Fig. 2.— Top: Integrated HCN map of NGC 253. Bottom: 
Integrated HCO + map of NGC 253. In both panels, the beam- 
size (3.9" X 2.7") is indicated in the bottom left corner. Con- 
tours are plotted starting at 5 Jybeam -1 kms -1 in steps of 5 
Jybeam -1 kms -1 . 



andlNguven-Q-Rieu et all ((1992) for both molecules and 
iPaglione et all (|1995l ) for HCN. 

The fluxes and luminosities given in Table [1] imply a 
HCN(1^0)/CO(1^0) flux ratio of 0.084 or a luminos- 
ity ratio of 0.14, and a HCO+(1^0)/CO(1->0) flux ra- 
tio of 0.077 or a luminosity r atio of 0.13, which is sim- 
ilar to other starbur sts (e.g.. iNguven-Rieu et all Il994t 
[ Jackson et"ai1 Il995t [Sorai et all l2002t iGao fc Solomon! 
l2004t iNakanishi et al.ll2005h . If only measuring the CO 
flux in those pixels where there is HCN emission, the CO 
flux decreases by 12%, indicating that there is relatively 
little diffuse CO emission present around the bright nu- 
cleus. 

For all three molecules we measure the peak sur- 
face brightness at the central position: CO(l— >0), 
2.4Jybeam" 1 ,; HCN(l-*0), 0.49 Jybeam" 1 ; and 
HCO+(1^0), 0.38 Jybeam" 1 . From this we determine 
the peak brightness temperature using the conversions 
6.95 JyK" 1 at 115.27GHz and 14.9 JyK" 1 at 89GHzfor 
our observations and derive the following temperatures: 
17K for CO(1->0), 7.3K for HCN(l->0), and 5.7K for 
HCO+(1^0). Thus, our CO data set a lower limit to 
the kinetic temperature > 20 K (Planck temperature). 
The CO(l— >0) brigh t ness t emperature is in agreement 
with IPaglione et al.1 ((2004), alt hough it is somewhat 
lower than the 36 K derived by iSakamoto et alj (2006) 
at comparable angular resolution. 

For several high redshift sources it has been found 
that lines from HCN appear to be narrower than the 
CO lines. This is the case for t he C l overleaf quasar 
(IWeifj et al.l 120031 : ISolomon et alj 120031: iRiechers et aT 



2006), APM 08279+5255 (fWagg et alJ 120051: 1 Weiss et al 
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Fig. 3.— left: APEX HCN(4^3) (thin black lines) and HCO+(4-»3) (thick grey lines) spectra maps towards the nucleus of NGC253. 
Spectra are shown on a velocity scale of -100 to 620 kms" 1 and on a T m t, brightness temperature scale from -0.1 to 0.4 K. Spectra have 
been observed on a regular raster with half beam spacing (9"). 
the central position. 



right: Enlargement of the HCN(4— >3) and HCO + (4— >3) spectra towards 



TABLE 1 
Resulting flux and luminosity 



Transition 


aperture 


/ 


V 








(arcsec) 


(Jy kms -1 ) 


(K kms -1 


pc 2 ) 


HCN(1^ 


0) 


20 


186 


9.4 X 


10 6 


HCN(1^ 


0) 


27 


208 


1.1 X 


10 7 


HCN(1^ 


0) 


Full map 


370 


1.8 x 


10 7 


HCN(4^ 


3) 


20 


1338 


4.5 x 


10 6 


HCN(4^ 


3) 


27 


1634 


5.5 x 


10 e 


HCN(4^ 


3) 


Full map 


2310 


7.7 x 


10 e 


HCO+(l 


-0) 


20 


178 


8.9 x 


10 fi 


HCO+(l 


-0) 


27 


203 


1.0 x 


10 7 


HCO+(l 


-0) 


Full map 


340 


1.7 x 


10 7 


HCO+(4 


-3) 


20 


1302 


4.1 x 


10 e 


HCO+(4 


-3) 


27 


1723 


5.4 x 


10 6 


HCO+(4 


-3) 


Full map 


2860 


8.9 x 


10 6 


CO(l-»Q) 


Full map 


4400 


1.3 x 


10 s 



Note. — The lines are detected at 20 — 30a. This is well 
within the 15% uncertainty of the flux calibration. 
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Fig. 4.— Global CO(l-tO), HCN(l^O) and HCO+(1->0) spec- 
tra of NGC253 derived from the OVRO data cubes. The CO spec- 
trum has been arbitrarily scaled down (by a factor 9.5). Each 
channel is 13.5 kms -1 wide. The same spectral region is shown as 
in Fig. 




2001 . and IRAS F1021 4+4724 (|Downes et all Il995t 
Vanden Bout et aI1l2004D . It has been argued that this 
indicates that the dense gas is more concentrated around 
the nucleus. Even though NGC253 does not match 
these high redshift galaxies in far-infrared luminosity, 
we extract the spectra for CO(l— >0), HCN(1— >0), and 
HCO+(1^0) (shown in Fig. gj). Here we find that the 
line widths are the same for all three lines. This agrees 
with the interferometric channel maps (Fig. [J), in that 
the emission emerges from the same cloud or system of 
clouds. 

For illustration purposes we also compare the distribu- 
tion of dense gas to Ha emission and 20 cm radio contin- 
uum in Figure[5] Strong Ha emission is seen in the dense 
gas regions traced by HCN(l->0) and HCO+(1^0). In 



the Ha image (from iLehnert fc Heckmar]|1995f ) the out- 
flow from the central r egion is very pr onounced. In the 
radio 20 cm map (from lUirostadl [2000h . we see that the 
radio emission traces the starburst region as well as the 
outflow. The 20 cm map compared to high-frequency ra- 
dio maps is particularly useful for comparison with the 
HCN and HCO + maps as the resolution is similar. The 
region towards the southwest is not seen in the Ha im- 
age, although radio continuum is detected at 20 cm, in- 
dicating that this likely is a highly obscured star forming 
region. 

From the line- free channels of the OVRO observations, 
we derive a 3 mm continuum image. The 20 cm emis- 
sion peaks at the center at a flux level of 0.1 Jybeam -1 , 
while the flux level in surrounding regions has dropped 
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Fig. 5. — Top: Ha image (greyscale) with HCN contours (contour 
levels as in Fig. [2}- Bottom: Radio continuum image at 20 cm, 
shown as greyscale with a logarithmic scale (again overplotted with 
the HCN contours). 
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Fig. 6.— HCN(1^0)/HCO+(1->0) ratio map of NGC 253. 
Values range from 0.7 (light grey) to 3 (black). The beamsize 
(3.9" X 2.7") is shown in the bottom left corner. 

to about O.OUy beam -1 . In the 3mm continuum map, 
the same is seen, where the flux level also peaks around 
0. Uy beam - 1 and also falls off by about an order of mag- 
nitude. The continuum emission as seen at 20 cm is dom- 
inated by synchrotron emission and the 3 mm emission is 
likely a combination of free-free emission and dust con- 
tinuum. A similar behaviour is seen in the Submillimc- 
tre Common-User Bolometer Array (SCUBA) maps from 
lAlton et al.l (|1999D , where the submillimeter observations 
trace the thermal dust emission. 

Figure [5] shows the HCN/HCO+ ratio map. Some of 
the pixels along the edges suffer from low signal-to- 
noise ratio, resulting in large uncertainties in the ra- 
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Fig. 7. — Histogram showing the distribution in the HCN(1— > 
0)/HCO+(1^0) flux ratio map (see Fig. |E}. The solid line his- 
togram represents the whole starburst region, while the dashed-linc 
histogram shows the distribution of the central region, where the 
continuum emission at 20 cm, 3 mm and 850 £tm peaks. The dotted 
lines show Gaussian fits to the distributions. 

tio in those pixels. We construct a histogram of the 
HCN(1^0)/HCO+(1^0) ratio map, shown in Figure [7] 
(solid line) , where the tails of the distribution are dom- 
inated by the low signal pixels. We perform a Gaus- 
sian fit to the histogram and find an average ratio 
(HCN(1^0)/HCO+(1^0)) = 1.1±0.3. To check for spa- 
tial variations of this value, we also show the histogram 
for the central region only (Fig. [Jj dashed line), where 
the continuum level is the highest. The ratio is 1.1 ± 0.2, 
while in the surrounding regions, where the continuum 
level has dropped by almost a factor 10, it is 1.1±0.4. We 
find no evidence for a change in the HCN/HCO + ratio 
across the regions with large variation in the continuum 
level (i.e., radiation field). 

4. DISCUSSION 

4.1. HCN and HCO + excitation 

A comparison of the 4-^3 to the 1— >0 line lumi- 
nosities of both molecules shows that the 4^3 tran- 
sitions are both subthermally excited independent of 



the aperture chosen (Lj 
U 



HCN(4^3) /-^HCN(l^O) 



= 0.45, 



hco+(4^3)/ L hco+(i^o) = °- 45 )- At a s P atial reso- 
lution of 20", the APEX data is directly comparable 
to previous HCN(4 — >3) measurements from the CSO 
(j Jackson et al.lfl995T ). Although the peak brightness tem- 
peratures of both measurements are roughly comparable, 
our integrated intensity is a factor of 2 lower than the 
measurements by Jackson et al (35 vs 71.4 K kms -1 ). 
A similar difference is apparent in the HCO + (4— >3) data 
(31 vs. 50.8 K kms -1 ). The APEX data in comparison 
to the Caltech Submillimeter Observatory (CSO) data 
have a much broader bandwidth and thus result in a 
much more reliable baseline subtraction and total fluxes. 
Given the quality of our new APEX data we adopt our 
values for the following excitation analysis. 

As the HCN(4^3) and HCO+(4^3) intensities are 
critical to constrain the excitation, we have reinvesti- 
gated the gas excitation by applying a series of spher- 



6 



Knudsen et al. 




2 3 4 5 6 7 
Rotational Quantum Number J . 




2 3 4 5 6 7 
Rotational Quantum Number J . 



Fig. 8. — HCO+ (right) and HCN (left) line SEDs towards the nucleus of NGC 253 in a 27" beam. The data points are from: 
INsniven-Q-Rieu et al.l Jl9§l, triangles'). iPaglione et~aT|[T997l and IPaglionel HMIt ) (HCN and HCO+, respectively, open square), and this 
work (Table 1; solid squares and circle). The lines in both diagrams are selected single-component LVG model fluxes for [n(H2), 7k; n ] 
combinations; solid line: 10 5,2 cm -3 , 50K; dashed line: 10 4 9 cm -3 , 150 K; dotted line: 10 5,5 cm -3 , 20K. The implied area filling factors 
are 1.4%, 2.0% and 1.6% for HCN and 0.9%, 1.6% and 1.1% for HCO+. 



ical, single-component, large velocity gr adient (LVG) 
model s. We use the collision rates from ISchoier et al.l 
(2005) and relative abundance ratios of [HCN/H2] = 
5 x 10- 9 an d pHCO+/H 2 1 = 1.6 x 10" 9 determined for 
NGC 253 bv lMartm etalf (|2006l) and a velocity gradient 
of 8 kms" 1 pc _1 . The latter is motivated from CO stud- 
ies which indicate that the a large velocity spread is re- 
quired in order to model the high 12 CO/ 13 CO line ratios 
obse rved towards NGC253 and other starburst galaxies 
(e-g JAalto et al.lll991t fPaglione et all 120011 : IWeifi et all 
2005). To compare the LVG predicted line brightness 
temperatures to the observations we convert the LVG 
line temperatures to integrated flux densities. Note that 
for fixed abundance ratios and a fixed velocity gradient 
the column density per velocity bin scales linearly with 
the density: N mo i/dV ~ X[mol] * n(H 2 ) / (dv / dr) . 

Figured] shows line spectral energy distributions 
(SEDs), i.e., rotational quantum number versus flux 
density, for HCN and HCO+ for the central 27" beam 
which allows us t o also in c lude the H CN(3— >2) and 
HCO+ (3-^2) from iPaglionl (fl997T l and IPaglione et all 
(|1997| ) to this plot. Both SEDs limit the possible range 



normalized HCN SEDs 
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in H 2 density to about 10 4 -10 s * 5 cm . Models with 
lower density fail to reproduce the observed 4^3 flux 
densities for both molecules. Similarly, higher densities 
over-predict the 3— >2 and 4— >3 lines even for very cold 
temperatures of 15 K. The kinetic temperature is poorly 
constrained from th e models (see also the LVG models in 
iJackson et al.lfl995l our Fig. [5], dotted and dashed lines); 
however, a kinetic temperature of 50 K is the most plausi- 
ble solution as it i s consistent with the dust temperature 
(|Melo et al.l I2002T) and also in agree ment with models 
of the CO SED (jGiisten et al.l l2006f ). We note that a 
tempe rature of ~ 100 K, as derived bv Bradford et all 
(2003), cannot be ruled out by the model. With this, 
a good match to both line SEDs in our fixed chem- 
istry framework is provided by an H2 density of 10 5 2 
cm~ 3 , a kinetic gas temperature of Tkin = 50 K, and an 
area filling factor of about 1% (Fig. [5] solid line). All 
solutions which match the data predict that HCN and 
HCO + lines are optically thick and subthermally excited, 



1 23456789 10 
Rotational Quantum Number J upper 

Fig. 9. — Comparison of the HCN line SEDs of selected lo- 
cal galaxies. The SEDs are shown for NGC2 53 (this paper; no 
data points plotted; Fig(SJ, M82 (soli d squa r es; [S caquist fc Fraverl 



2000), Arp220 (open squares; Greve et al. (2009) and references 
therein), NGC6240 (open triangles: IGreve et a!Tl|2006l ) and ref- 
erences therein [For the SED fit to NG C6240 we ha ve used the 
HCN(l^O) flux measurement offS olomon ct al. (1992) as the more 
recent value by[Grcvc ct al. (2006) does not provide a good fit in 
conjunction with the HCN ( 3— >2) and HC N (4— >3 ) measurements.]), 
Mrk231 (soli d triangles; Papadopoulosl l|2007l' ), and NGC4945 
(open circles; IWang et al.l 120041^ The HCN line SEDs are nor- 
malized by their HCN(1— >0) flux density. 



even in the ground transition of HCN and HCO + . For 
the 50K model we get an excitation temperature T ex of 
26K for HCN(l^O) and of 42K for HCO+(1-^0) and all 
higher transitions have lower T ex . For a very cold solu- 
tion of 20K {n H2 = 10 5 - 5 cm^ 3 ) the HCO+(1->0) emis- 
sion would be close to being thermalized. Although the 
area filling factors for HCN and HCO + differ slightly for 
a given n# 2 - 2\i n LVG input, the good agreement of 
the model predicted SEDs for HCN and HCO+ with the 
observations suggests that the emission from HCN and 
HCO + arise from a common volume. 
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Fig. 10.— Flux ratios of HCN(1^0)/HCO+ (1^0) plotted as 
function of 20 cm flux for each pixel. The three crosses embedded 
in a diamond are taken from the central position very close to the 
peak of 20 cm emission, where the angular HCO+ flux distribution 
hints at a divergence from the angular HCN flux distribution. 

We repeat the analysis for the central position of 
NGC253 after convolving the data to a resolution of 
20". For our LVG model with n(H 2 ) = 1Q 5 2 cm" 3 and 
Tkin = 50 K, the results are consistent with those derived 
from 27" resolution. The only difference is the area fill- 
ing factor which increases for higher resolution. If we 
apply the same LVG model to the total flux within the 
maps covered by the OVRO and APEX observations, the 
LVG predicted HCN and HCO+ (4-v3) fluxes are slightly 
higher than the observed values which suggests that the 
HCN and the HCO + excitation is somewhat lower in the 
outer most regions of the data cubes. At a fixed temper- 
ature of Tki n = 50 K, however, a small decrease of the H 2 
density to n(H2) = 10 51 cm -3 is sufficient to bring the 
model prediction and the observations into agreement. 

In summary, the constant HCN/HCO+ flux ratio 
across the starburst disk would require a highly unlikely 
combination between varying abundances and excitation 
conditions. Hence it is safe to conclude that there is no 
significant abundance gradient across the disk. 

4.2. Comparison to excitation conditions in other 
galaxies 

We have used our new HCN observations as well as 
data from the literature to compare the excitation of 
the dense gas via the HCN SEDs in NGC253, M82, 
NGC4945, NGC6240, Mrk23I and Arp220 (references 
are given in the caption of Fig. [9|). We restrict our- 
selves here to HCN because it provides the most complete 
data set among tracers of the dense gas for these active 
galaxies. All galaxies except Mrk23I have been mea- 
sured in HCN(l->0), HCN(3->2), and HCN(4^3). For 
Mrk23I, only the HCN(I^O) and HCN(4^3) have been 
published in the literature. Figure [9] shows a comparison 
of the HCN SEDs of the six galaxies normalized by their 
HCN(l^O) flux density. This comparison shows that 
among these sources NGC 253 exhibits one of the high- 
est HCN excitations. Given that the HCN emission in 
NGC253 is subthermally excited, this comparison shows 
that subthermal HCN is a common property of star form- 
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Fig. 11. — HCN/HCO + luminosity ratio vs. the far-infrared lu- 
minosity for spiral/s t arbursts galaxies ( Nguvcn-Q-Ric u et aljl992t 
llmanishi et' al. 2004, 2006; Wang ct al. 2004J) and LIRGs/ULIRGs 
IIGracia-Carpio et ajj|200ejl . We have included our measurements 
for NGC 253, the recent high redshift detection o f HCO+ in the 
Cloverleaf quasar (Ricchers ct al. 2006 1 HCN fro m Solom on et al.l 
I2003T) and in the APM0827 9+5255 lIGarci'a-Burillo et al.l 120061 
HCN from [Wagg ct al. 20Q5|), and indicate the ratio for t he cir- 
cumnuclear disk of the Milky Way (Christopher ct al. 2005). Note 
that for all galaxies the 1 — » transition was used, expect for 
APM08279+5255 where the 5 — > 4 transition was used. For the 
high-z objects, a standard concordance cosmology is used, with 
H = 73kra s- 1 Mpc -1 , Q M = 0.24, and U A = 0.72 flSpergel et al.l 
[20031. 120071 ). 

ing galaxies. The only source where the HC N excitation 
exceeds that of NGC253 i s Arp220 (see also iGreve et al.l 
120061 lPapadopouiosll2007n . We not e however, t h at the 
HCN line intensities of Arp220 in IGreve et all (|2006l ) 
(which have been used to fit the HCN SED of Arp220 
in Fig. [9]) are higher by a factor of ~ 1.5 — 3 compared to 
other studies (ISolomon et al.lll992t IWiedner et al.1 [20021 
ICernicharo et alj 12006) . Therefore, the available data is 
not conclusive and does not allow to rule out that HCN 
is subthermally excited also in Arp220. Interestingly, 
NGC6240 and Mrk231 show a lower excitation of the 
dense gas phase compared to NGC253 (Fig. [9]) despite 
their much higher FIR luminosities (Fig. [TT|) . 

4.3. Possible interpretation of £hcn/^hco+ var "l a ^ ons 
between galaxies 

While our results show that the HCN/HCO + ratio is 
statistically constant across the starburst disk, we do in- 
spect the details. In Figure [10] we plot the HCN/HCO+ 
flux ratios as function of the 20cm fluxes for each pixel. 
Towards higher 20cm flux, there are three data points 
which hint at a higher HCN/HCO + ratio (this increase 
is due to a relative decrease of the HCO + emission). 
They correspond to the peak in the 20cm flux. This 
is also towards the peak of the hard X-ray emi ssion as 
measured with Chandra ([Strickland et al.ll2000t ), where 
it is possible that there is an absorbed active galactic 
nucleus (AGN) pre sent in the very center of NGC253 
(|Weaver et al.| [2002') . While these data points are within 
the scatter, we speculate that the small increase in the 
HCN/HCO+ ratio could be a consequence of the pres- 
ence of an AGN b ased on the relations suggested by 
iKohno et al.l (j2001) . This, however, can only be resolved 
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Fig. 12.- LVG predicted i' HCN{1 ^ 0) / £' HGO+ (1 _ 0) 
function of the HCN (top) and the HCO+ abundance (bottom). 
The lines show LVG models with a H2 density of 10 4,5 cm -3 (light 
grey squares), 10 5 cm -3 (black circles) and 10 55 cm -3 (grey tri- 
angles). The kinetic temperature has been fixed to T^ n = 50 K, 
the adopted velocity gradient is 8 kms - 1 pc _1 for all models . 

when future instruments with larger resolving power and 
higher sensitivity will become available. 

The line ratio between HCN and HCO + has been stud- 
ied for several classes of star for ming galaxies, mostly 
through single-dish observa t ions (|Nguven-Q-Rieu et all 
1992t llmanishi et all 12004 teOOGU Gracia-Carpi oet all 



2006), where ratios have been found in the range 0.5- 



2.7. iGracia-Carpio et al.l (|2006[ ) presented a data set of 
HCN/HCO+ for luminous IR galaxies (LIRGs; L m > 
10 11 L Q ) and ultra-luminous IR galaxies (ULIRGs; 
Lm > 10 12 Lq). They present a possible correlation 
between the HCN/HCO + ratio and the IR luminosity. 
If, however, other measurements from the literature, in- 
cluding those of less luminous star forming galaxies and a 
recent high-redshift measu r ement (iNguven-Q-Rieu et al.l 
[1991 llmanishi et aT][2006l l200i iNakanishi et all 120051 ) 
are included in such a plot, there remains little evi- 
dence for such a cor relation with Ltr (Fig- HU see also 
iRiechers et"alll2006| ). 



The large number of galaxies with ^hcn/-^hco+ 



> 1 



seems surprising at the first glance since the emission in 
both lines is likely optically thick in all galaxies (see also 
iGreve et al.l 1200 6). and the critical densities of HCO + 
are lower than those for HCN. From this one would 
expect that the HCO + (l— >0) line is brighter than the 
HCN(l^O) line if both molecules trace the same gas. 
However, as both molecules are subthermally excited, 
their line intensity depends on the underlying chem- 
istry despite their high optical depth. This is exem- 
plified in Figure [I2l where we show the LVG predicted 
^hcn/-^hco+ ra ti° as a function of reasonable HCN and 
HCO + abundance for different H2 densities. From the 
plot it can be seen that a decrease of the HCO + abun- 
dance can indeed explain even the highest observed line 
ratios independent on the underlying density of the gas. 
Alternatively, an increase of the HCN abundance will 
also increase the -^hcn/-^hco+ ra ti°- The latter process, 
however, is only effective for H 2 gas densities below a few 
times 10 5 cm~ 3 , as for higher densities HCN(1— >0) will 
become thermalized and, therefore, insensitive to abun- 



dance variations. Given the gas density derived by us 
for NGC253 (10 52 cm -3 ), both processes can effectively 
modify the observed line ratio. 

A decrease of the HCO + abundance due to the ioniz- 
ing field produced by co smic rays has been proposed by 
iSeaquist fe Fraverl (|2000f ) . Cosmic rays ionize H 2 , leading 
to the production of H;j~, which reacts with CO to form 
HCO + . The abundance of HCO + is thus affected by ra- 
tio of cosmic-ray ionization rate and gas density. While 
a higher ionizing flux favours the production of HCO + , it 
also increases the number of free electrons, which leads to 
a higher probability for recombination. At a gas density 
of 3 x 10 4 cm~ 3 , an ionizing field comparable to that of 
the Galaxy may already be strong enough to significantly 
decrease the abundance of HCO + due to dissoc iative re- 
combination of H^" (e.g. Phi llips fc Lazio|1995f). 

T heoretical model l ing b y iMeiierink fc Spaani (|2005l ) 
and iMeijerink et aTJ (|2007f ) suggest that for modest col- 
umn densities Nh < 10 cm -2 the presence of an 
embedded X-ray source, such as an AGN, can enhance 
the HCN/HCO+ and HCN/CO abundance ratios rela- 
tive to a pure UV source, such as young massive stars 
in a starburst. Note, however, that for larger column 
densities the reverse would be true. There are hints of 
enha nced HCN/HCO+ ratios in nearby Seyfert galaxies 
(e.g.. lKohno et aLll2001[ ). Nonetheless, similar enhance- 
ments a£einot_Jound fo^_the far- infrared bright Clover leaf 
QSO (^lechers et al.l 12001 . 

An alternative explanation for the observed variations 
of the ^hcn/^hco+ ra, ti° could be that HCO + and 
HCN are not arising form the same volume. From 
an analysi s of the dense gas properties in Arp220 and 
NGC6240, IGreve et al.l (|2006f ) recently concluded, that 
the HCO + emission is dominantly emitted from regions 
with 10 times lower density than HCN. In this picture 
the ^hcn/^hco+ ra ti° would then be explained by the 
relative area filling factors of the denser cores and the 
less dense surounding medium in combination with the 
excitation conditions in both gas phases. This picture, 
however, is not supported by our observation of NGC253, 
where the excitation of both molecules can be explained 
by a single density-temperature combination. 

5. CONCLUSIONS 

We have presented maps of HCN(l-^O) and 
HCO + (1^0) obtained with OVRO and observations of 
HCN(4^3) and HCO+(4^3) obtained with APEX of 
the starburst disk of NGC253. We find that the spatial 
distribution of the two molecules is very similar, even 
when compared on a channel- by-channel basis. We 
find the ratio between HCN(l^O) and HCO+(1^0) is 
1.1 ±0.2 across the whole circumnuclear starburst region 
independent of the strength of the central radiation field. 
The emission from the HCN(l-O) and HCO+(1-0) 
transitions, both indicators of the dense molecular gas, 
trace regions that are non-distinguishable within the 
uncertainties of our observations. 

From an excitation analysis we find that both 
molecules are subthermally excited and that their emis- 
sion lines are optically thick. The constant line ratio im- 
plies that there are no strong abundance gradients across 
the starburst disk. A consequence of both molecules be- 
ing subthermally excited is that the line intensity is sen- 
sitive to the underlying chemistry despite their high op- 
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tical depth. This finding may also explain the variations 
in £hcn/-^hco+ between different star forming galaxies. 

The advent of the Atacama Large Millimeter Array 
(ALMA), Expanded Very Large Array (EVLA), and 
Square Kilometer Array (SKA) will open up improved 
prospects for observations of dense molecular gas phase 
in galaxies through molecular line transitions such as 
HCN and HCO+ at low and high redshifts. As both 
HCN and HCO + appear to trace the same dense gas re- 
gions where the formation of the stars are taking place, 
observations of these two molecules will be important 
probes in future projects aimed at studying galaxies dur- 
ing phases of the build-up of their stellar mass. 
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